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Abstract 

The Causal Dynamical Triangulation model of quantum gravity (CDT) has a 
transfer matrix, relating spatial geometries at adjacent (discrete lattice) times. 
The transfer matrix uniquely determines the theory. We show that the measure- 
ments of the scale factor of the (CDT) universe are well described by an effective 
transfer matrix where the matrix elements are labeled only by the scale factor. 
Using computer simulations we determine the effective transfer matrix elements 
and show how they relate to an effective minisuperspace action at all scales. 



1 Introduction 



Minisuperspace models of the universe provide us with simple quantum mechan- 
ical models of fluctuations of the scale factor of the universe. In the simplest 
models one assumes spatial homogeneity and isotropy. Classically this implies 
that we can write 

ds 2 = -dt 2 + a 2 (t)dn 2 (1) 

where dfl 2 is the line element of a homogeneous and isotropic three-dimensional 
space. Assuming space is compact it is S 3 . Under these assumptions the dynam- 
ical variable is the scale factor a(t) of the universe and the quantum field theory 
of the universe is reduced to quantum mechanics of a single variable a(t). 

The 4D Causal Dynamical Triangulation model (CDT) of quantum gravity is 
by construction a (regularized) quantum field theory model which is compatible 
with spatial homogeneity and isotropy (for reviews see P, E]). It uses the path 
integral formulation and assumes there exists a foliation in (proper) time. When 
the average over all geometries of this kind is performed one indeed finds that 
the average geometry, i.e. what one naively would think is closest to a "classical 
geometry", can be described by a line element of type (jl]). Even more, it turns 
out that to a good approximation the scale factor a(t), measured as the average 
of the third root of the spatial volume at (proper) time t is well described by the 
simplest minisuperspace action. Contrary to the standard use of minisuperspace, 
where one postulates the reduction to geometries described by a metric like Q 
and then only quantizes the single degree of freedom a(t), the CDT discussion 
of a(t) is exact (to the extent that CDT describes quantum gravity). The a(t) 
entering in the CDT discussion is obtained by including all geometries in the 
path integral. It is natural to ask the following question: how well does the 
simplest minisuperspace action describe the CDT data generated by Monte Carlo 
simulations. The tool we will use when trying to answer this question is the 
transfer matrix. 

The time foliation present in CDT provides us with a transfer matrix. The 
geometries considered in the (regularized) path integral are piecewise linear ge- 
ometries constructed in such a way that at discretized times t n the spatial slices 
are triangulations of S 3 . The transfer matrix relates a given spatial (piecewise 
linear) geometry at time t n to a given spatial geometry at time t n +i- In ordinary 
(Euclidean) lattice field theory reflection positiveness of the transfer matrix en- 
sures a unitary time evolution. The CDT transfer matrix has a similar property 
[3J [2] [j] In this article we will analyze an "effective" transfer matrix of CDT to 

1 The idea of a time foliation and the fact that there exists a related unitary time evolution 
are features which CDT shares with Hofava-Lifshitz gravity [3]. However, no spatial higher 
derivative terms are explicitly added to the action, like in Hofava-Lifshitz gravity, and it is 
possible that fixed points of the lattice theory can be identified with the non-trivial UV fixed 
points conjectured in the asymptotic safety scenario suggested by Weinberg [5] and investigated 
in HE]. 



2 



be defined below. 



2 CDT and the CDT transfer matrix 



The use of piecewise linear geometries was introduced in the context of general 
relativity by Regge [9j as a natural tool to work with a discretized version of 
the Hilbert-Einstein action, but without the use of coordinates. The idea was to 
approximate a given smooth geometry by a continuous, piecewise linear geometry. 
This piecewise linear geometry is uniquely defined by a triangulation where the 
length of the links are given. The observation by Regge was that the standard 
Einstein-Hilbert action in D dimensions, 



for such piecewise linear geometries has a geometric interpretation as a sum over 
deficit angles of the D — 2-dimensional sub-simplices in the triangulation. While 
the idea of Regge was to approach a given classically smooth geometry by a 
sequence of suitable piecewise linear triangulations, a different use of piecewise 
linear geometries was made in the formalism of Dynamical Triangulations (DT) 
|lUj . Although we have no mathematical rigorous definition of the path integral 
over geometries, it is natural, in analogy with the path integral in ordinary quan- 
tum mechanics, to assume that the summation over geometries will involve not 
only smooth geometries but all continuous geometries. A subclass of these ge- 
ometries is the piecewise linear geometries and a further subclass is the piecewise 
linear geometries defined by triangulations obtained by gluing together equilat- 
eral D-simplices such that they form a manifold of fixed topology. In DT the 
assumption is that this set of geometries is in a suitable sense dense in the set 
of continuous geometries when we take the link length a to zero, and in this way 
the link length a will act as an ultraviolet cutoff, just like in ordinary lattice 
field theory. Further, the natural signature of space-time in the DT formulation 
is Euclidean. We will assume this is the case in the rest of this article. That 
this approach works in principle, i.e. that it is able to reproduce a continuum 
quantum field theory which is diffeomorphism invariant, is well documented for 
D = 2. Two-dimensional Euclidean quantum gravity coupled to a conformal field 
theories with c < 1 can be solved analytically both in the continuum [Til EE] and 
using the formalism of DT [T3] and agreement is found. 

In higher dimensions DT was studied numerically, using Monte Carlo sim- 
ulations both in three dimensions [U] and four dimensions [Jjj]. However, no 
convincing continuum limit has been obtained so far in higher dimensions [16J, 
and this was one of the motivations for changing the class of triangulations used 
in the path integral in CDTrl In the CDT formalism one sums over geometries 




2 It has recently been suggested that there might be a continuum limit of DT which belongs 
to the same universality class as the CDT theory described below 
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with a (proper) time foliation. In principle one starts out with space-times with 
a Lorentzian signature (contrary to the situation in DT) and the foliation is 
in proper (Lorentzian) time. However, each piecewise linear geometry used in 
the CDT path integral allows a rotation to Euclidean proper time. The set of 
Euclidean geometries we obtain in this way is a subset of the DT Euclidean ge- 
ometries and this restriction seemingly cures some of the higher dimensional DT 
diseases, while in two dimensions the relation between the restricted theory and 
the full DT theory has been worked out in detail: one obtains the CDT theory 
from the DT theory by integrating out all baby universes (which results in a 
non-analytic mapping between the coupling constants of the two theories), and 
(somewhat surprisingly) one can restore the DT theory from the CDT theory by 
the inverse mapping [T?l ITS] . Using four-simplices (which is the case having our 
attention in this article) as building blocks one can, for suitable choices of bare 
coupling constants, observe a four-dimensional (Euclidean) universe [IH]. For 
these choices of coupling constants the shape of the universe is consistent with an 
interpretation as an (Euclidean) de Sitter space, at least as long as one looks at 
the scale factor [20J . This is the region of coupling constants which will have our 
interest (for other choices of the coupling constants one obtains more degenerate 
configurations [2~T]). 

An interesting feature of the CDT model, rotated to Euclidean signature, is 
that it possesses a transfer matrix (3J, [2] . At (discrete) time t n we have a spa- 
tial hypersurface with a spatial geometry characterized by a three-dimensional 
triangulation T 3 (t n ). At time t n+ i we have spatial geometry defined by an- 
other three-dimensional triangulation Ts(t n+ x). We assume for simplicity that 
the topologjj^] of the spatial triangulations is that of S 3 . In CDT we sum over 
all four-dimensional triangulations of the "slab" between t n and t n+ \ compatible 
with the topology S 3 x [0, 1] and such that each four-simplex which "fills" the 
slab has subsimplices which are (sub)simplices of T 3 {t n ) as well as T 3 (t n+1 ). This 
leads to 4 types of four-simplices in the slab: type (4, 1), (3,2), (2,3) and (1,4), 
where the numbers denote the number of vertices in T 3 (t n ) and T 3 (t n+ i) respec- 
tively. The number N^'^^tn) of (4, 1) simplices in the slab is equal to the number 
N 3 (t n ) of three-simplices in T 3 (t n ), and similarly the number of (1,4) simplices 
in the slab is equal to the number N 3 (t n+ i) of three-simplices in T 3 (t n+1 ). The 
transfer matrix A4, i.e. the amplitude between the T 3 (t n ) and T 3 (t n+ i), is now 
given as the sum over all such triangulation, 



where the summation is over all four-dimensional triangulations of a slab, with 
boundary triangulations T 3 (t n ) and T 3 (t n+ i), CV 4 is the order of the automorphism 

3 S 3 is chosen for simplicity. We could have chosen any spatial topology. The important point 
in the assumption is that the topology is not allowed to change from one spatial hypersurface 
to the next. 




(3) 
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group of the triangulation T 4 , and where S[Ti] is the Regge action of the four- 
dimensional triangulation of the slab. 

The transfer matrix is defined on the vector space spanned by the set 7% of 
three-dimensional triangulations. This space is infinite dimensional and has the 
natural scalar product 

(T\V} = ^ r 6 T , T ,, ^\T)C T (T\=I, T,T'e% (4) 

T T 

where Ct is the order of the automorphism group of the triangulation T. 

The transition amplitude for a three-dimensional triangulation T to develop 
into a three-dimensional triangulation T' after t tot + 1 (integer) time steps is 

(T'\M^ +1 \T) = ^(T'|>l|T ttot )C Tttot (T iiot |>l|T iiot _ 1 )---C ri (T 1 |M|T). (5) 

We will define the partition function corresponding to t to t time steps as path 
integral with periodic boundary conditions after t to t time-steps: 

Z ttot = trM*** = J2Ct(T\M^\T). (6) 

T 

This is the partition function we have used in our computer simulations. The 
measurements performed so far, using Monte Carlo simulations, have been con- 
centrated on the measurement of the scale factor, or more conveniently the three- 
volume n ti = Ns(ti) at the spatial slice at time U, as well as the correlation 
between the three- volumes at time ti and time tj. These "observables" can be 
expressed using the transfer matrix Ai. The probability of measuring the spatial 
volume n ti at time ti is given by 

p'"°'(n k ) = -!- Yl orvw~m = ^g^ - m 

In ([7j %(n ti ) denotes the subset of three-dimensional triangulations where the 
number of three-simplices is n ti and p{n ti ) the projection operator on the subspace 
spanned by these triangulations: 

P(n)= £ \T)C T (Tl p{nf = p{n). (8) 
Similarly the correlator between n tl and n t2 , separated by At = t2 — t\ is given 

by 

P^(n tl ,n t2 ) = -L Ct^M^-^T^Ct^M^T,) 

ttot Tl er 3 (n tl ) 

T 2 e%(n t2 ) 

tr pjntJM^-^pMM^ 

tiM*** ' { ' 
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and this expression can clearly be generalized to multi-correlators. 

As mentioned in the Introduction it was possible to explain the observed 
distributions P ttot (n tl ) and P ttot (n tl ,n t2 ) using a simple minisuperspace model. 
How does the concept of a minisuperspace labeled by states n ti , i = 1, . . . ,t tot 
relate to the transfer matrix A4 which is defined on the much larger space spanned 
by the vectors T? Let us define the "effective" transfer matrix M by 

\n)(n\M\m)(m\ = ^ \T n )C Tn (T n \M\T m )C Tm (T m \. (10) 

T n eT 3 (n) 
T m eT 3 (m) 

(n\M\m) represents the average of the matrix elements (T n \Ai\T m ) like 

(n\M\m) = (M) n , n := T X r E (T n \M\T m ), (11) 

T n £Ti(n) 
T m eT 3 (m) 

where M n denotes the cardinality of T n . Af n grows exponentially with n. 



In (10) and (11) it is misleading to think of the "state" \n) as the (suitably) 
normalized sum of the Af n vectors \T n ) (although by doing so one would of course 
obtain the correct expectation value (M) n ,m using such vectors). Such a vector 
would again be a single vector located in the M n - dimensional space spanned by 
the |T n )'s. It is more appropriate to think of the "state" associated with n as 
arising from a uniform probability distribution of states \T n ) and in this way to 
think of p(n) as the associated density operator. However, once we have reduced 
our consideration to the matrix (n\M\m) we are of course free to find eigenvectors 
for this matrix and expand them in the abstract basis \n), and we will indeed do 
that. 

The statement that we can use the matrix (n\M\m) as an effective trans- 
fer matrix is the statement that the standard deviation of the Af n Af m numbers 
(T n \A4\T m ) is sufficiently small. In fact the difference between trAi 2 and trM 2 
can exactly be expressed as a sum over deviations squared for each n, m: 

trM 2 -trM 2 = ^ ]T (VC^C^(T n \M\T m ) - (M) n , m y (12) 

n,m T„e73(n) 
T m eT 3 (m) 

In the following we will assume that we can work with an effective transfer 



matrix (n\M\m). Eq. (10) is an attempt to define this effective transfer matrix 
from first principles and in principle one can check by computer simulations if 
it is a good approximation. We will here take the pragmatic attitude to assume 
there exists such an effective transfer matrix and use it to analyze the computer 
generated data. The consistency of this analysis is indirectly evidence that an 
object like (n\M\m) provides a good approximation of our data. Thus we will 
use the "effective" version of (J^J) - (Joj) : 

p(n) = \n)(n\, (13) 
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t , \ tr p(n t )M ttot , . 

where p(n) should be distinguished from p(n). 

In particular we can measure the matrix elements (n\M\m) up to a normal- 
ization by considering t tot = 2. We have: 

P m M= MMMMMM (16) 

V V trM 2 V J 

This method requires a major change in our general computer program which 
assumes t tot > 3. We updated it but we were not completely convinced that the 
new version is stable although it gave exactly the same results as the method we 
finally used. For t tot = 3,4 we get: 

tr M 6 

(4) = MWMMWM (18) 

V 1> A) trM 4 V ) 

From the measurements of P^ (711,712) and P^ 4 \ni^n^) we can determine the 
matrix elements (n\M\m) up to a normalization: 

(n\M\m) = C -j= ^== \ (19) 

There is nothing magic about the above choice. One could have chosen t tot = 4 
and tt t — 6 and formed the combinations 

p(4)Kn2) = m 
pm{nuni) = MMMV!> (21) 



from which one can again extract (n\M\m) like in ([19]). We have indeed checked 
that measurements of P^(n 1 ,n 2 ) and Pw (711,714) lead to the same M matrix as 
extracted from measurements of P( 3 )(n 1; n 2 ) and P^{ni^ n 3 ), up to a normaliza- 
tion. 

In earlier work we have shown [201 122] that the following minisuperspace 
actionlll 

err n 1 (n t +i-n t ) 2 fn t + n t+l \ 1/3 n t + n t+ i 



(22) 



1/3 

4 In fact we used slightly different form of the potential terms: fi ( nt+ ™ t+1 j — A " t-f '"' — 
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1/3 

/m t — An t . This parametrization was more convenient to extract the parameters of the action 
from the measured covariance matrix of volume fluctuations. In this article we implement a 
modified form ( 22 ) which better fits our data. 
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describes well the measured (n t ) and the fluctuations (n t n t /) — (n t ) (n t >) in the bulk 
where n t is large. The effective action (|22|) suggests that the effective transfer 



matrix 



(n\M\m) oc e 



(n-m) i / n + m \ 
n + m "T^y 2 ) 



1/3 



\ n+m 



(23) 



is a good approximation in the bulk. We will in the following try to determine 
the transfer matrix from the data, also in the range where n t is not necessarily 



large and we will try to improve the expression (23) 



3 How to perform the computer simulations 

The simplest version of the discretized CDT theory has three parameters, two 
related to the cosmological constant and the gravitational constant, and an ad- 
ditional parameter which controls the asymmetry between the edge lengths in 
the spatial and time directions. This latter parameter seems not to be a genuine 
coupling constant since it just labels the different length assignment of spatial 
and time-like links. The action used is still the Einstein-Hilbert action (as for- 
mulated by Regge for piecewise linear geometries), adjusted for this asymmetry. 
However, because we study the theory in a truly non-perturbative region of cou- 
pling constant space the effective action is determined by a competition between 
the classical action used and a contribution coming from the measure term. The 
contribution from the measure term is "entropic" in nature: it counts the num- 
ber of configurations with the same action and is thus independent of the other 
parameters. Effectively this promotes the asymmetry parameter to a genuine 
coupling constant (see j2] for a detailed discussion). 

In the numerical simulations the topology of the manifold is assumed to be 
S 3 x S 1 with periodic boundary conditions in the (Euclidean) time, as mentioned 
above. The four-simplices used to construct the simplicial manifolds of CDT are 
characterized by their position in spatial and time directions. As also mentioned 
above we have four types of four-simplices: (4, l)-simplices, with four vertices at 
time t and one vertex at t+ 1, (3, 2)-simplices with three simplices at t and two at 
t + 1 and the "time-reversed" (1, 4)-simplices and (2, 3)-simplices. All simplices 
of a particular type are identical. 

The discretized (Regge) Einstein-Hilbert action becomes extremely simple 
because we are essentially only using the two kinds of building blocks to construct 
the four-dimensional triangulation T [3J [2j : 

Sr[T] = -(«o + 6A)iV + /€ 4 (iV (4,1) + iV (3 ' 2) ) + A (2A^ (4 - 1) + iV (3 ' 2) ) (24) 

where N is the total number of vertices in the triangulation, jV^ 4,1 ) the total num- 
ber of type (4,1) plus (1,4) simplices and jV^ 3 ' 2 ) the total number of simplices 
of type (3,2) plus (2,3). Kq, and A are the (bare) dimensionless coupling 



8 



constants obtained by the discretization of the continuous action ([2]). k is pro- 
portional to the inverse bare gravitational constant, k 4 related to the cosmological 
constant while A is related to the asymmetry between the spatial and time-like 
links. A = corresponds to spatial and time-like links having the same length. 
An additional geometric parameter is the length ttot of the periodic time axis. 
The partition function 

Z = e ~ SR[T] (25) 
TeT 4 

has a critical value Kf lt (Ko, A), depending on k q and A, such that Z is divergent 
for k 4 < K,f lt . The existence of this critical value is reflecting the fact that 
the number of triangulations with a fixed number of four-simplices iV 4 grows 
exponentially with iV 4 . In principle we want to fine tune k 4 to this critical value 
since we really want a limit where N± — > 00. In practice the simulations have so 
far been carried out by keeping N4 (or iV^ 4 ' 1 )) fixed. In this way we have been 
trading the coupling constant k 4 with iV 4 and the partition function Z(N±) is 
related to Z(k 4 ) by a Laplace transformation: 

Z(k 4 ) = Y,e~ KiN4 Z(N 4 ). (26) 

The phase diagram now depends on kq and A and we refer to [2J [21] for a detailed 
discussion. Here we will be working in the interesting, so-called de Sitter phase 
where we, for a given (large) iV 4 and sufficient large t tot , observe a (Euclidean) de 
Sitter universe, i.e. a four-sphere where the temporal extension is proportional to 
N^ 4 while the rest of the time extension (assuming t to t is large enough compared 
to iV 4 1/4 ) is a stalk of almost no spatial extension. Presumably this stalk only exists 
because our computer algorithm does not allow the spatial extension to shrink to 
zero. In Fig. [llwe have shown a typical situation with t tot = 80, iV^ 4 ' 1 ) = 160000 
and we observe a bulk region (the "blob" , approximately from t = 20 to t = 60) 
where nt, the three- volume, i.e. the number of tetrahedra at the time-slice t, is 
large, and the rest is the stalk region where nt is very small. Fig. [T] shows both 
the average over many configurations and a typical configuration which appears 
in the path integral. When taking the average over many configuration we align 
the center of mass of the blobs (see [20] for a detailed discussion). 

We can also measure the probability distribution P ttot (rit) of n t in the blob 
for a given t in Fig. [I] It is shown in Fig. [2] (left figure). It is well approximated 
by a Gaussian distribution around the mean value (nt). This is in contrast to the 
situation in the stalk where the probability distribution splits in three families 
[23J, as shown on the right part of Fig. [2j 

Figs, [l] and [2] are based on computer simulations of the type mentioned above: 
jy(4,i) j s kept fixed. Technically this has been done by adding a term e(A^( 4,1 ^ — 
j\r( 4 ' 1 )) 2 to the action, e being a suitably small parameter: 

Sr ->• S R + e(N^ - N^) 2 . (27) 
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Figure 1: The distribution n t of the three-volume, i.e. the number of (4, 1)- 
simplices, when we are in the de Sitter phase. The data is obtained for k = 2.2, 
A = 0.6. A single (typical) configuration as it appears in the path integral is 
shown by the blue line, and the average distribution is indicated by the red line. 

This term ensures that iV^ 4,1 ) is going to fluctuate not too far from iV^ 4 ' 1 ). The 
precise value (A^ 4,1 )) depends on the choice of k±. We now want to study the 
transfer matrix. However, the structure of the transfer matrix is incompatible 
with a global constraint of this type, so we have to change the updating procedure. 
We have done this in two different ways. The first way is to drop the constraint 



term e(A r ^ 4 ' 1 ' ) — iV^ 4 ' 1 )) 2 and only use the discretized Einstein-Hilbert action (24). 
The way to obtain an average (N^ 4 ^) is to fine tune n 4 to K^ lt '. The closer k 4 
is to the critical value the larger (N^' 1 ^). In practice this fine tuning can be 
difficult and the larger the system the more difficult the fine tuning. Thus we can 
and will only use it for small systems. For larger systems we apply a different 
strategy which also constrains the value of N^ ,l \ but which is compatible with 
the transfer matrix structure: we change the global constraint imposed on iV^ 4 ' 1 ) 



in (27) to a local constraint in t: 

ttot 



S R ^S R = S R + s J2( n t - n vol f. (28) 



t=i 



Of course this constraint will drastically change the profile (nt), since n t will 
now fluctuate around n vo i. Thus we will have different transfer matrices Ai and 
M, and different probability distributions P{n%, n,2, ■ ■ ■)■ The new probability 
distribution for n t is shown for various n vo i in Fig. [3j 

However, we can reconstruct some of the probability distributions associated 
with the action Sr if we know it (i.e. measure it) for the action Sr. The proba- 
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Figure 2: Left figure: Probability distribution of the volume n t in the blob (t = 29 
in Fig. IT]), Right figure: Probability distribution in the stalk (1 < t < 17 in Fig. 

0. 



bility for measuring (ni, n 2 , . . . , n ttot ) is given by 

(n 1 |M|n 2 )(n 2 |M|n 3 ) • ■ ■ (n tt jM\m) 



P^(n 1; n 2 ,...,n ttot ) 



tr M ft tot 



and is directly related to the distribution without the volume fixing term, 



P ttat {n u n 2 ,...,n tiot )^P ttot {n u n 2 ,...,n tt Je e(ni n ^ )2 • • • e 
We calculate the transfer matrix M in the same way as M: 

p( 3 )(ni = n,n 2 = m) 



<nt tot -n vol ) 2 



(n\M\m) 



P( 4 )(nx = n,7i3 = m) 



(29) 



(30) 



(31) 



To calculate the original transfer matrix M we have to cancel the volume fixing 



term, which is easily done: From equations (|28h , (29) and (30) we obtain 



(n\M\m) 



e \e(n-n vol ) 2 



(n\M\m)e 



\e(m-n vol ) 2 



(32) 



For each choice of n vo i we observe n t with some approximate Gaussian distri- 
bution centered around n vo i, where the width depends on our choice of e, and we 
use the associated probabilities to construct (n\M\m), as described above. To re- 
construct the matrix in a larger region of the n-space we have to merge data from 
different n vo i regions. Since the matrix is determined only up to a normalization, 
the way to do this is to make sure there are regions of overlap between the n t 
distributions and in these regions choose a suitable calibration procedure such 
that we can merge the data. We will later describe how this is explicitly done. 
Needless to say the procedure we are employing here is a kind of multi-canonical 
Monte Carlo method (see [21] for a review). 
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Figure 3: Probability distribution of the volume n t , for various n vo i = 
200,400,600,1000,1400,1800 and 2200 (from left to right). For all ranges, the 
simulations were performed with K4 = 0.322 and e = 0.00002. 

4 The effective action at large three-volumes 

We can measure the transfer matrix (n\M\m) for large n,m as described above. 
As already noted it is well approximated by the matrix M^: 



(n\M {th) \m) =Afe- L ^f {n ' m \ 
where the effective Lagrangian is 



J eff 



[n, m) 



1 
f 



(n — m) 2 
n + m — 2ri( 



+ /i 



n + m 



1/3 



A 



n + m 



(33) 



(34) 



We now ask how well? We make a best fit of the parameters T, /x, A and M (fixing 
uq = 0). The measured M, M'" 1 ' from (33) as well as their difference, are shown 
in Fig. [4] for n vo i = 1400 and for the n t range 1200 < n t < 1600. 

The values of parameters T, \i and A for different values of n vo[ , obtained from 
the best fits of the matrix (33) to the measured matrix M are presented in 



Table [TJ Again no is chosen to be zero. 
4.1 The kinetic term 

To get a better estimation of the parameters associated with the effective action 



(34), we first try to fit only to the parameters of the kinetic term which is by far 
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Figure 4: The left figure: The transfer matrix M for range 1200 < n t < 
1600. The matrix is optically indistinguishable from the fitted theoretical transfer 
matrix 

jlf (t/0. The right fi gure: The difference between M and disappears 
in the numerical noise. The measurements were performed for k,q = 2.2, A = 
0.6, k 4 = 0.3220 and e = 0.00002. 





n t range 


r 


n 


A 


600 


400 - 820 


25.71 ±0.01 


18 ±1 


0.05 ±0.01 


1000 


780 - 1220 


26.00 ±0.01 


17 ± 1 


0.05 ±0.01 


1400 


1180 - 1630 


26.10 ±0.01 


13 ±1 


0.04 ±0.01 


1800 


1580 - 2040 


26.08 ±0.01 


26 ± 1 


0.07 ±0.01 


2200 


1980 - 2440 


26.05 ±0.01 


19 ±2 


0.05 ±0.01 



Table 1: The values of and A for different n vo i, 
M^ th > to the measured M. 



obtained from best fits of 



the dominating term from a numerical point of view. We do that by keeping the 
sum of the entries, i.e. n + m, fixed such that the potential term is not changing. 
In this way we can try to determine T and even n which we had put to zero 
in the fits mentioned above in order not to have too many fit-parameters. The 
matrix elements for constant n + m = c show the expected Gaussian dependence 
on n (see Fig. [5J: 



(n\M\m) = (n\M\c — n) = M{c) exp 



(2n - cf 



T • c- 2n n 



(35) 



where the terms in the effective action which only depend on c are included in 
the normalization. 

We expect the denominator of the kinetic term k(c) to behave like k(n + m) = 
r • {n + m — 2no). As shown on Fig. [6] this is indeed true and the parameter T is 
common for all ranges. Fig. [6] presents measured coefficients k(c) for various c's 
and ranges of n t denoted by distinct colors together with a linear fit. The best 
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Figure 5: (n\M\c — n) plotted as a function of n for c 
(dots). Gaussian fits are drawn with a line. 



2000, 2800 and 3600 



linear fit gives Y = 26.07 and no 
for separate ranges of nt- 



-3, which is consistent with results obtained 



4.2 The potential term 

The potential part of the effective Lagrangian may be extracted from the diagonal 
elements of the transfer matrix 



1 



L e ff(n,n) = - log(n\M\n) + c(n vol ) = - (fm 1/3 - Xn) 



(36) 



However, because of different normalizations of the transfer matrices for different 



n,n) 



ranges (hence the dependence of the constant c{n vo i) on n vo i), the fit of £ e // 
to the transfer matrix data cannot be performed in a straightforward way. The 
transfer matrices have first to be merged properly via a scaling procedure, i.e. by 



adjusting the c constant in (36). This is done in the following way. For example, 
for n vo i = 1400 the range of spatial volumes for which we measured the transfer 
matrix is n t = 1180 . . . 1630, while for n vo \ = 1800 the range is n t = 1580 . . . 2040. 
Thus, there is a non-vanishing intersection n t = 1580 . . . 1630 for which elements 
of both matrices were measured. We scale the second matrix, so that the mean 
value of the diagonal elements on the intersecting region is equal for both matrices. 
After applying this procedure for successive ranges, we finally get scaled transfer 
matrices which can be merged. The result of such merging is shown on Fig. [7j 
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1000 2000 3000 4000 

c 



Figure 6: The coefficient k(c) in the kinetic term as a function of c = n + m 
(different colors denote different ranges), and a linear fit k(n + m) = T ■ (n + m — 
2no) (red line). 



which shows the diagonal elements of the scaled transfer matrices together with 
a fit of form (36), where we took V = 26.1. 



4.3 A global effective action fit 

Summing up, the effective action determined via the transfer matrices is strikingly 



well described by eq. (|34|). As in the last subsection we can merge the scaled 

pression (34) to the 
7, fi = 15.0 and A = 0.046. 



matrices for all the choices of n vo i (see Fig. [8]), and fit expression (34) to the 
aggregated data. The best fit gives T = 26.17, n 



Method 


r 


n 




A 


Cross-diagonals 


26.07 ±0.02 


-3± 1 






Diagonal 


(26.07) 




16.5 ±0.2 


0.049 ± 0.001 


Full fit 


26.17 ±0.01 


7±1 


15.0 ±0.1 


0.046 ±0.001 


Previous method* 


23 ± 1 




13.9 ±0.7 


0.027 ±0.003 



Table 2: The values of T,no,fj, and A fitted in different ways. *We also present 
the parameters of the effective action extracted from the covariance matrix of 
volume fluctuations in our earlier work 
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1000 1500 

n 



2500 



Figure 7: log(n|M|n) of the scaled transfer matrix (dots, different colors denote 
different ranges) compared with the fit of the potential term —L e ff (red line, 
which stops at n — 400). 



We summarize the results of fitting the parameters of L e ff in different ways in 
Table [2] As long as we are concentrating on the large (bulk) values of n t the 



various data clearly do not allow us to improve the expression (34). 

For comparison we also present the parameters of the effective action measured 
indirectly from the covariance matrix of volume fluctuations. This method is 
based on the analysis of the effective propagator around the semi-classical solution 
|23j . The parameters of the action from the previous method agree quite well with 
those measured directly from the transfer matrix. The small difference may result 
from slightly different parametrization of the potential term in the effective action 
(see footnote 4). 

4.4 Miscellaneous 

The spectral decomposition of matrices presents us with an interesting way to 
compare the measured transfer matrix (n\M\m) with the "theoretical" matrix 



(n\M^\m) given by eq. (33) with parameters obtained from a best fit, as de- 
scribed above (in this case the fit from the n vo \ = 1400 data). We obtain very 
good agreement between eigenvalues and eigenvectors for the two matrices. The 
first 6 eigenvalues are presented in Fig. M and the first 4 eigenvectors are shown 



in Fig. 10 



As a final check of the consistency of the effective transfer matrix with data we 
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Figure 8: Logarithm of merged and scaled transfer matrices measured for different 



measure the probability distribution of a spatial slice volume P( ttot \n) in terms 
of the transfer matrix M given by 



P 



ttot 



n) 



(n\M ttot \n) 
tr M**°* 



(37) 



We use M instead M, because the power M ttot involves summation over all pos- 
sible volumes n which are not accessible until we suppress them with the term 



e € ( n Uv ° 1 ^ (which is present in M). We compare the theoretical expectation (37) 
with the measured P ttot (n) for t tot = 3 and t to t — 4. The comparison is shown on 



Fig. 11 and the error is of order 0.02%. 



5 The transfer matrix for small three-volumes 



We have seen that the effective transfer matrix is very well described by the 
simplest effective Lagrangian (34). However, it is natural to expect that the 



corresponding effective action is only a first approximation. Appealing to isotropy 
and homogeneity one would expect that a potential R k (t) term, where R(t) refers 
to the scalar curvature of the three-dimensional space at time t, translates into a 
term ( nf+ ^ t+1 ) 1 2k ^ 3 _ The leading term R(t) is already part of the effective action 
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Figure 9: The first 6 eigenvalues of the measured transfer matrix (n\M\m) cal- 
culated for n vo i = 1400 (dots) and the similar eigenvalues for the fitted transfer 
matrix (n\M^\m) (line). 

where it appears as the term ( nt+ ^ tf+1 ) 1 ^ 3 . For the data coming from large n t we 
have seen that the difference between our measured M and coming from 

1 /3 

the effective action with only the ( Ti ' t+ ^ t+1 ) term is already at the noise level. 
From these large n t data we have no chance with the present statistics to study 
higher k corrections. Thus we now turn to the small n t region. The difficulty of 
analyzing the small n t region is that we might encounter discretization effects as 
is apparent from Fig. [2j The study of this region started in [23] , but the present 
approach offers the great advantage that we can perform high statistics study of 
small systems, while in the earlier studies the interesting region was a small part 
of a larger system and thus the statistics becomes less good. 

We want to measure the effective transfer matrix by measuring P' tot (nj, rij) 
for t tot = 3,4 as we have already done, but now for small n t , i.e. for much 
smaller systems. We also want to check as well as possible that the concept of 
an "effective" transfer matrix actually works. While it seemed to work well for 
large n t , it is not obvious that this will remain true for small n t : physics might be 
different and discretization effects might also spoil such a picture. Our systems 
will be so small that we do not have to introduce the auxiliary transfer matrix 
M and the volume fixing parameter n vo i. Thus we only use the Regge Einstein- 



Hilbert action Sr, eq. (24), to generate the probability distributions P ( - 3 \ni,n 2 



and P^{ni,ns) and construct (n\M\m) from eq. (19). 
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Figure 10: The first 4 eigenvectors of the measured transfer matrix (n\M\m) (cal- 
culated for n vo [ = 1400 (dots)), and the first 4 eigenvectors of the corresponding 
theoretical transfer matrix {n\M^ th '\rn) (lines). 

5.1 Eigenvectors analysis. 

Given the transfer matrix M we can perform a spectral decomposition in terms 
of eigenvalues A, and (orthonormal) eigenvectors \a.i): 

M = J2 X M)(oa\. (38) 

i 

Since the measured M is only determined up to a normalization, we will assume 
Ai = 1 and |Ai| > |A 2 | > .... 

If the gaps are significant between the first, second and third eigenvalues , it 
is clear that the large ttot limit of Z{t to t) and the large At limit of the simplest 
correlation functions P ttot (n t ,n t +At) will be completely dominated by the first 
two eigenstates: Thus, in the limit where t to t ^ At ^> 1 we have (recalling that 
Ai is normalized to 1) 

Z{Uot) = ^) Uot ~ 1 ( 39 ) 

i 

P ttot (n) = -^Yl ^ ( n \ a i) 2 ~ H«i) 2 ( 40 ) 

i 

P u °*(n t ,m t+At ) = |(^(A,) A '(n|« l )(« l |m))(^(A,)' tot - A *(n|a i )(a i |m)) 

i i 

~ (n\a\} 2 (ai\m) 2 + A Ai (n|a2)(«2|^)(^|ai)(ai|^) 
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Figure 11: The probability distribution P ttot (n) (dots) and ^f^j^ (line) for 
Uot = 3 (left) and t tot = 4 (right). 

The average can be written as: 

(n)t tot « ^n(n|ai) 2 (41) 

n 

For the correlator we have: 

(n t m t+At ) - (n t )(m t+ At) ~ A^* nm(n|o!2) (a 2 1 m) (n|ai) (oti |m) . (42) 

n,m 

and the long distance behavior is an exponential fall off e _AtAt , fi = — logA2/Ai 
(where we have reintroduced Ai for clarity). 



i.o - 

0.8 - 
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k 4 = 0.33 



k 4 = 0.3223 



0.0 
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Figure 12: Eigenvalues of the transfer matrix M for different values of K4 ap- 
proaching K c p l Rj 0.3222. Note that according to our normalization the biggest 
eigenvalue is set to one. 

How well are these approximate relations satisfied? First we observe that 



there is indeed a clear gap between the first eigenvalues, as illustrated in Fig. 12 
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The closer k 4 is fine-tuned to the smaller the gap, but even very close to 
the critical value we observe a clear gap. That we have gaps even at nf lt just 
illustrates that we indeed consider small systems. 

Next we ask how large ttot has to be in order that the approximations made in 
eq. (40)-(42) are valid. That can of course be read off from the eigenvalues and 
already for t tot = 4 the approximation is very good for (n) ttot . For the correlator 
one has of course to consider larger t tot . in Fig- 13 we have shown for t tot = 12 the 
expected exponential decay with exponent log(A2/Ai) compared to the actually 
measured correlator. The agreement is very good even for small At where it is 
not obvious that ignoring the eigenvectors with eigenvalues smaller than A2 is a 
valid approximation. 



100 



Empirical, T = 12 
Empirical full CDT 
Theoretical 



10 



■- 
a 

o 



0.1 




Figure 13: The "theoretical" correlator (42) of spatial volumes between time- 
slices separated by At and calculated including the first two eigenvectors of the 
transfer matrix M calculated for k 4 = 0.3223 (green line). The correlator decays 
exponentially with At (log scale). There is a very good agreement with the 
correlator measured directly in simulations for T=12 (red) and in "full CDT" 
stalk range (blue). The bars indicate measurement errors. 



5.2 The "Full-CDT" approximation. 



We have also compared our approximate large t tot probability distributions P(n) 
(eq. 40) with the data taken from the stalk range of full CDT simulations (in- 
cluding the blob and the stalk range). The distributions approach very well "full 



CDT" measurements as K4 tends to critical value (Fig. 14). The validity of the 



transfer matrix model is further confirmed by the behaviour of the correlator 
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(n t m t +At) — ( n t) ( m t+At) measured directly in the stalk range of "full CDT". As 
illustrated in Fig. 13 the measured correlator falls off as e~^ At . The parameter fi 
is well explained by the ratio of the first two eigenvalues of the transfer matrix 
M calculated for k 4 closest to the critical value. 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




Full CDT 

K4 = 0.3223 
ka = 0.35 
ka = 0.40 




Figure 14: Theoretical 1-point probability distributions calculated with the 
'largest' eigenvector of the transfer matrix M approach empirical probability 
measured in "full CDT" (red line) as n± tends to critical value. Theoretical dis- 
tributions were computed for k 4 = 0.40, 0.35 and 0.3223. For k 4 = 0.3223 (blue 
circles) the agreement is very good. 



6 The effective action for small three-volumes 

In principle the matrix M presented above would allow us to determine an effec- 
tive Lagrangian for small n t , i.e. even in the stalk range of the CDT configurations, 
precisely as we did for large n t : 

S eff = J2Leff(n t ,n t+1 ), (n\M\m)=Afe- L ^f^ m \ (43) 
t 

However, we are confronted with the existence of three families of states, as is 
apparent in Fig. [2J We can however define a reduced matrix M performing a 
summation over the three families, i.e. 

M = UMU T (44) 
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where the rectangular matrix U has a form: 



U 



1 1 1 
1 1 1 
1 1 1 



The elements of the matrix M behave much more smoothly and can be analyzed, 
using the effective action idea. We normalize M by choosing its largest eigenvalue 
to be one. 

Let us assume that L e jj has the form: 



L eff (n,m) 



(n — m)" 



+ v(n + m) 



k{n + m) 

where the functions k(-) and v (•) are to be determined. 



(45) 



<n|M|c-n> 

0.025 r 
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0.000 




100 150 200 250 300 350 



Figure 15: (n\M\c — n) plotted as a function of n for c = 222 (blue), c = 282 



(red), c = 342 (green). Gaussian fits (46) are presented as solid lines 



We follow the same strategy as when n t was large: first we analyze matrix 
elements for constant n + m = c in order to keep the potential term v(n + m) 



constant. One observes a Gaussian dependence on n (Fig. 15): 

(2n - cf 



(n\M\m) = (n\M\c — n)— J\fexp 



k(c) 



v(c) 



(46) 



Fitting (|46|) for different c's one can easily check that the kinetic coefficient 

(47) 



k(c) is linear (Fig. 16). Therefore one can write: 

(n — m) 2 



J eff 



[n, m) 



1 

r 



+ v(n + m) 



n + m — 2tiq 

where: vQ = TvQ. The best fit of T and no is presented in Table [3 



23 




Figure 16: Kinetic coefficients k(c) (blue) and a linear fit (red). 



In the above we recognize the familiar kinetic term present in the effective 
action for the blob range. Let us also test the assumption that the potential part 
is similar by analyzing the diagonal elements of M where the kinetic term is zero. 



From (43) and (47) one obtains 



log(ra|M|n) 



T(2n) + log M 



(4* 



log<n|M|n> 




n 



100 200 300 400 



Figure 17: Logarithm of the diagonal (n\M\n) (blue dots) compared with the fit 



of the potential term (49) (red line). 



Let us assume 



v(2n) = —Xn + fin 1 / 3 + 5n 



(49) 



inspired by the bulk expressions already derived and where we have included a 
new term Sn~ p , inspired by earlier remarks about effective powers of the three- 
dimensional curvature R(t). We present the fit of the logarithm of the diagonal 
elements (n\M\n) in Fig. 17 and the parameters of the fit in Table [3 
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1 (XL cLlilc Lcl 




JJlUU 


1 


27.2 ± U.l 


OCT op o 

25. r — 2b. 2 


n 


± 1 


— o h f 


n 


34 ±2 


13-30 


A 


0.12 ±0.02 


0.04 - 0.07 


5 


(4 ± 7) x 10 4 




P 


3± 1 





Table 3: Fitted parameters of the effective action for the stalk (50). For compar- 
ison we also present estimates of parameters of the effective action for the blob 
calculated from the large n t simulations (see Table [I] and [2]) . 



It is quite surprising that the same effective action is still present in the stalk, 
despite the volume behavior seems, at the first sight, quite different from that in 
the blob range. It is even more surprising that the parameters of the fit agree 
quite well with the effective Newton constant V measured in the blob range. 
The parameter p is slightly bigger but of the same order of magnitude as the 
potential coefficient from the blob range. Only A which is related to the size of 
the dynamically created universe is quite different, but that should be no surprise 
since A semiclassically is related to the size of the universe. Finally the value of 
p is difficult to explain from the point of view of higher powers of R(t) which 
should give p = 2k/3 — 1, k = 2, . . ., but also it should be mentioned that it is 
not very well determined from the fits. 

Summing up: 



ostalk 
^eff ■ 



[n t - n t+1/ 



n t + n t+1 -2n 



+P 



nt+riu 



1/3 



-A 



n t +n t+1 



h5 



n t +n t+1 



(50) 

Having determined by a best fit the parameters of the effective action we can 
calculate the "theoretical" transfer matrix M^ th "> using (43) and (pOh. To appre- 



ciate the quality of this approximation we present a plot of six lowest eigenvalues 
of the measured M and the "theoretical" matrix M^ th ' as well as the comparison 
of their six lowest eigenvectors. In each case the continuous line corresponds to 



the AfW (Fig. 18-19). 



7 Discussion and conclusions. 

CDT comes with a transfer matrix (T\Ai\T'). The way CDT is defined allows 
us to measure certain distributions, say P ttot {n t ) and P ttot (n t ,m t+ At), of three- 
volumes n t . These distributions have an exact definition in terms of the transfer 
matrix M. and the density matrices p(n) which are projections onto the subspace 
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Figure 18: First six eigenvalues of the measured M (blue dots) and the theoretical 
matrix (red lines). was calculated using the effective action for the 

stalk range ((501). 




Figure 19: The first six eigenvectors of the measured M (blue dots) and "theo- 
retical" matrix (red line). M (t/l) was calculated using the effective action 



for the stalk range (50). 



of three-dimensional combinatorial triangulations T(3) of S 3 spanned by the tri- 
angulations with n tetrahedra. (see eqs. ([7])- ([9])). While the transfer matrix is 
defined on the large vector space spanned by the elements in T(3), the actual data 
coming from Monte Carlo simulations seem to allow for a much simpler descrip- 
tion in terms of an "effective" transfer matrix M, only labeled by abstract vectors 
\n) referring only to the three- volume. Not only that: basically over the whole 
range of n t the data are described by a transfer matrix which can be represented 
as 



(n\M\m) = Af e - L *ff( n > m \ 



(51) 
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where L e ff(n,m) is given by 



j eff (n,m) 



r 



(n — m) 2 ( n+m\ 1 ^ 3 , / n+m\ „ / n+m^ 



n+m — 2n \ 2 

"(52) 

and with a corresponding effective action 



The last term in (52) unfortunately is not very well determined. For large n t 
we can not really observe it. The first terms seem to fit that data perfectly 
with the present statistics. For small n t one can detect a term like 5 (*4p) p , 
but as mentioned, still p is not well determined, and in addition the value we 
obtain depends on the specific merging of the three different distributions one 
observes for small n t . Thus we cannot really claim that we have a result which is 
discretization independent. We are caught in an unfortunate dilemma: we want 
to go to small n t in order to observe this term, which indicates corrections to the 
simplest minisuperspace action. However, taking n t small also brings us into the 
region where discretization effects are likely to be important. It is possible that 
one can find a window where n t is small enough for the term to be observed via 
high statistic measurements, but where n t is large enough to avoid discretization 
effects, but we have not yet pursued this in a systematic way. 

This discussion highlights an important advantage of the present method: 
since t tot is small, we are effectively simulating much smaller systems than in the 
traditionally "full" CDT computer simulations. In this way we can actually obtain 
measurements of high statistics with relatively moderate computer resources and 
in a finite amount of time. 

The amazing accuracy with which the effective transfer matrix seems to be 



described by eq. (51) indicates that one obtains a good approximation to the 



partition function by writing: 

Z ttot = tr M ttot = e~ Se!![nt >\ (54) 



where S e ff[n ti ] is the effective action (52)- (53). We have strictly speaking only 



shown that this expression is a good approximation for some special values of 



the bare coupling constant of the Einstein-Hilbert action Sr, given by eq. (24). 
However, without much doubt any choice of the bare coupling constants which 
places us well inside the so-called de Sitter phase will allow for a description in 
terms of an S e ff[n ti ], just with different T, \i and A. Let us assume that this is 
also true in the two other phases which have been observed in the "full" CDT 



theory. If this is the case one can actually use S e ff[n tj ] and the expression (54) 



to study the phase structure of CDT. This is of course much easier than using 
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the full system. This is precisely what has been done in a recent paper 
Seemingly one obtains a good qualitative description of the CDT phase diagram 
(and also new interesting phase structures), corroborating the conjecture that the 
functional form (52)- (53) might be sufficient to describe CDT for all choices of the 
bare coupling constants of the Einstein-Hilbert action (24). Checking this is an 
obvious task for the future. However, an even more interesting application of the 
multi-canonical Monte Carlo simulation method developed here is that it might 
allow us to investigate the CDT phase transitions in more detail. A possible UV 
scaling limit of the CDT theory has to be associated with these phase transitions. 
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